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We study the appearance of resonantly coupled optical modes, optical necklaces, in Anderson 
localized one-dimensional random superlattices through numerical calculations of the accumulated 
phase. The evolution of the optimal necklace order m* shows a gradual shift towards higher orders 
with increasing the sample size. We derive an empirical formula that predicts m* and discuss the 
situation when in a sample length L the number of degenerate in energy resonances exceeds the 
optimal one. We show how the extra resonances are pushed out to the miniband edges of the 
necklace, thus reducing the order of the latter by multiples of two. 



Wave interference phenomena play a crucial role in 
transport properties of various physical systems from 
periodic [1] to disordered @, Among these, origi- 
nally studied for electronic systems, Anderson localiza- 
tion seems the more intriguing one [4] . It predicts a phase 
transition form metallic-like conductivity to an insulat- 
ing regime, when the transport can come to halt with 
increasing the randomness. On the other hand, in a one- 
dimensional (ID) disordered system, such a transition 
happens when the sample size L extends beyond the so- 
called localization length £. In such conditions, the wave- 
functions become localized within an extension £ and de- 
cay exponentially with distance I. Being essentially an in- 
terference phenomenon, Anderson localization has been 
studied as for electromagnetic and acoustic waves [H , as 
well as for degenerate atomic gases 0|. Very recently, 
Anderson localization of optical waves in the microwave 
regime has been demonstrated in experiments on ID ran- 
dom multilayer dielectric stacks 0, [§[ • 

Initially, it has been widely accepted that the conduc- 
tivity (transmittivity) of a disordered chain is mainly 
supported by states which are closer situated to the sam- 
ple center [9|. Later, this was questioned, since for long 
enough specimen, the states in the center possess signif- 
icantly reduced probability to support a two-step hop- 
ping transport through these states. In late 80's, Pendry 
[10[ and Tartakovskii et al. [ljj, independently, argued 
that the conductivity in such systems should be domi- 
nated by so-called necklaces - few homogeneously dis- 
tributed through the sample states, degenerate in en- 
ergy, and coupled resonantly to delocalize and extend 
through the chain. The number m of resonant states 
forming a necklace, was calculated to scale as L 1 / 2 , while 
the probability of their occurrence was shown to drop 
as exp(— L x l 2 \ thus predicting them to be increasingly 
improbable events in long samples [10]. Therefore, for a 
certain sample length, a trade-off between the expected 
number m and their occurrence probability would deter- 
mine the optimal order of the necklace. 

In this Letter we study the appearance of optical neck- 
laces in finite ID random superlattices through numeri- 
cal calculations of the accumulated phase and follow the 



evolution of the optimal necklace order when increasing 
the sample size. We suggest an empirical formula which 
predicts the optimal necklace order m* and discuss the 
situation when in a sample length L the number of de- 
generate in energy resonances exceeds m* . We show that 
the extra resonances, which are spaced by a distance less 
than the optimal one for a certain sample length, are cou- 
pled strongly enough to split and be pushed out to the 
miniband edges of the necklace, thus reducing the order 
of the latter by multiples of two. 

We studied random binary multilayer stacks composed 
by AT layers of A and B-type dielectric materials (re- 
fractive indices chosen to be tia = 1.3 and n# = 2.1). 
Positionally random sequences were generated by giving 
equal weights to each type of layer. Two types of ran- 
dom multilayer stacks were investigated. In one type, the 
physical thicknesses oIa and oIb of the constituent layers 
were generated randomly through the sequence, thus a 
complete randomness was achieved. For the other type, 
oIa and ds were chosen to be resonant at some wave- 
length Ao (quarter- wave layers), therefore, the condition 
rtA&A — tlboIb = Ao/4 was considered. The quarter- 
wave condition makes the system not perfectly random 
and, in particular, maximizes the constructive interfer- 
ence effects at the resonant wavelength. 

In general, in a ID dielectric multilayer stack, the 
phase suffers a tt shift each time it meets a resonance. 
Therefore, the phase measurement is a valid tool to iso- 
late the spectrum singularities in either periodic [12[ or 
random systems |8[ . This way, in Ref. [8[ , localized states 
and, more importantly, optical necklaces have been iden- 
tified through interferometric measurements on Anderson 
localized ID random superlattices. When m > 2 degen- 
erate in energy states couple weakly to form a necklace 
[13[, the accumulated phase increases smoothly through 
the transmission band built up by these latter, summing 
up to mir. On the contrary, when strong coupling oc- 
curs, the phase suffers single tt jumps through well de- 
fined spectrally separate peaks. 

Hence, through this study, the order of the neck- 
laces has been assigned to the multiples of tt in the 
smoothly changing phase shift through the transmittance 
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FIG. 1: (Color online) The logarithm of integrated trans- 
mission versus sample length for two cases of non-resonant 
(O) and resonant (■) layers. The solid and dashed lines are 
linear fits to the data for non-resonant and resonant layers 
cases, correspondingly. The fits give localization lengths of 
^8.5 //in and ~ 10.7 /im respectively. The histograms on the 
right show the counts of the order of observed necklaces for 
five different thicknesses of samples with resonant A/4 layers. 



band. Firstly, random structures, with non-resonant 
layer thicknesses, were studied and their transmission 
spectra were calculated through a standard transfer ma- 
trix method [l4|. In Fig.[T]the logarithm of the integrated 
transmission over a wide spectral range (1 -i-2 /im), aver- 
aged over 200 realizations per point, is plotted versus the 
sample thickness (open circles) [15j. We observe, firstly, 
that < InT > decreases as L -1 , as it is expected in the 
diffusive regime [2|, where the states are still extended. 
Over about 30 jam a smooth transition from extended 
to localized regime occurs and the sample transmission 
drops linearly with further increase in the sample thick- 
ness. The localization length of £ ~ 8.5 /im is obtained 
from the slope of < InT >= —j. 

In order to study the properties of optical necklaces, 
one has to go through many realizations of disorder, be- 
cause of their nature to be extremely rare events. For 
this, we decided to look for necklaces in binary multilayer 
stacks with resonant quarter- wavelength layers. We were 
motivated by the fact that, due to the resonant layer con- 
dition, these provide much higher probability to possess 
many resonances at Ao wavelength. 

As it is seen from Fig. [TJ the transmission of such real- 
ized samples (full squares) shows a clear deviation from 
the one of fully random system. Namely, the higher-on- 
average transmission in this case is supported by much 
frequent occurrences of optical necklaces at Ao and be- 
comes more pronounced for thicker samples. We have 
verified the order of the necklaces in a small frequency 
window around the resonant wavelength for five different 
sample lengths. For this, 300 realizations per each cho- 
sen thickness have been studied and the necklaces have 
been identified through the calculated phase. The right 
panels in Fig. [1] report the counts of different order neck- 



laces. We observe that at smaller thicknesses (panel A, 
L « 103 /im) the major number of necklaces is of second 
order and fewer third order ones occur. With increasing 
sample size, the number of second order necklaces dimin- 
ishes and the third order ones firstly start to dominate at 
L « 184 /im, obtaining a maximum at L « 220 /im. With 
a further increase in sample size (panel E, L « 253 /im) 
higher orders shoulder starts to rise in the necklace counts 
histogram, indicating that the samples become enough 
thick to fit more than three degenerate resonances. 

The order of the necklace, which fits optimally inside 
the sample length L to have the highest and most com- 
pact transmission band, can be empirically calculated as 

m* = + 1, (1) 

where |_ J stand for the /foor-function and denote the 
greatest integer part of the expression inside. Here, the 
parameter a counts the distance (in the units of £) be- 
tween two neighboring resonances efficiently coupling to 
form a necklace, while /3 considers the coupling of first 
and last resonances to the environment (see Fig. [2]). In 
reality, the step function of Eq. [T] considers the order 
of the most frequent necklace from the histogram of all 
observed orders. In our case, a reasonable fit to the nu- 
merical results was found with a = 7 and /3 — 2.2. 

The power law scaling of the order of necklaces in 
disordered systems and the very low probability of the 
higher orders occurrence predicts a trade-off scenario, 
which should reduce the number of actually observed 
ones [10]. In the following, we address this issue to under- 
stand how this reduction happens and reveal interesting 
physical picture of the phenomenon. 

Let us consider the case of an optimal rath order 
necklace stretching through sample length L. The most 
compact and, in the meanwhile, high- transmission band 
formed by this should occur when ra resonances, degen- 
erate in energy, are distributed homogeneously through 
the sample to favor similar coupling between neighbor- 
ing states. This situation is an ideal one and is very im- 
probable for higher order necklaces. Since the resonances 
occur randomly in the depth of the sample, it becomes 
more probable that some of them will appear at a dis- 
tance closer than the optimal separation I* = a£ between 
two neighboring cavities. When this happens, the spa- 
tially close resonances couple strongly and a larger mode 

Ipi; r=ct£, 
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FIG. 2: A simplified model of formation of an optimal-order 
necklace through the sample length L. 
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FIG. 3: (Color online) The light intensity distribution inside 
a periodic superlattice of four coupled cavities (left panels) 
and the corresponding transmission and phase spectra (right 
panels) . All the spectra are plotted in a the narrow frequency 
range around z/o = c/Ao = 200 THz. (a) homogeneously dis- 
tributed cavities (an artificial ra=4 necklace) form the high- 
est transmission band through which the phase accumulates 
smoothly a 4tt increase, (b)-(d) moving the two central cavi- 
ties closer to each other, the light transmission around Ao at- 
tenuates by orders of magnitude and the smooth phase jump 
counts only 2tv, reducing the necklace order. 



repulsion pushes them farther from the miniband center. 

We examine the dynamics of such a situation on an 
example of a periodic photonic structure, where for sim- 
plicity four resonant cavities are coupled through Bragg 
sequences in such a way that a high transmission band 
of an artificially built necklace is formed. The choice of 
a periodic structure allows us to manipulate the spatial 
positions of various cavities, which is an impossible task 
for random systems. Figure [3] shows the light intensity 
distribution through the photonic structure (left panels) 
and the corresponding transmission spectrum together 
with the phase (right panels). In Fig. [3]A. the ideal sit- 
uation is presented, when the resonances are distributed 
homogeneously through the sample. A high transmission 
miniband is formed (left panel, black line) through which 
the phase (red) suffers a smooth increase by 47r. 

Next, we move spatially the two central cavities closer 
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FIG. 4: (Color online) An example of a necklace formed in a 
218 /xm-thick random system. The intensity distribution map 
(left) shows six bright clumps, while the phase (right) counts 
four 7r's, indicating the formation of a fourth order necklace. 



to each other. The intensity distribution plot in Fig. [3f> 
shows the effect of stronger coupling between these two 
resonances: two well defined resonances (symmetric and 
anti-symmetric modes) rise out of the miniband core 
due to the mode repulsion and the channel intensity de- 
creases. The transmission of the miniband peak drops by 
four orders of magnitude, and the smooth change in the 
phase now sums up to 2tt. If we move these cavities 
closer (Fig. [3]0,D), the peak transmission drops down 
to ~ 10 -8 . Note that when excluding the two central 
cavities (while maintaining the same sample length) the 
transmission drops further down by almost three orders 
of magnitude to have a peak transmission of ~ 10 -11 . 

The intensity map in Fig. [3p contains an interest- 
ing feature, to which we would like to draw a special 
attention. We observe three maxima instead of two 
through the original transmission miniband, while the 
phase counts always a 2tt jump. We stress, that the cen- 
tral intensity maximum is not due to a resonance present 
at Ao, but is built by the overlap integral of the trans- 
mission lines of the symmetric and anti- symmetric states. 
Therefore, it is correct to consider the necklace at Ao to 
be of second rather than of third order one. 

In order to catch such a situation in a random system, 
we have performed a number of realizations of a 950 layer 
sample (218 /im-thick). In Fig. [H we show such an event. 
In the right panel six well-resolved intensity maxima are 
observed, while the smooth change of the phase sums up 
to 47r. The intensity distribution map, and, therefore, 
the electric field profile at the miniband center frequency 
shows many bright and other less intense clumps, but 
the real number of resonances, directly involved in the 
light transport through the miniband is less and can be 
exactly counted through the phase variation. 

Thus, when the number of randomly built degenerate 
cavities exceeds the optimal order of the necklace for a 
certain sample length, the extra cavities couple strongly 
to be pushed out to the edges of the miniband, reducing 



4 



the order of the necklace state and its peak transmission. 
We stress that this reduction is always a multiple of two, 
independently on the fact whether even or odd number 
of extra cavities are coupled strongly. This is because a 
strong coupling of even or odd number of cavities results 
in a dip or a peak in the spectral line at the resonant 
wavelength. Therefore, at Ao, only multiples of 2ir are 
filtered out from the phase jump. Nevertheless, the ex- 
tra cavities play an important role in the formation of 
necklaces, since their non-zero overlap reduces the reflec- 
tivity of inter-cavity regions and contributes positively to 
link the resonant tunneling transport through the neck- 
lace miniband. 

Along this, in a fully random system, another inter- 
esting situation can happen. Suppose, that next to a 
necklace band centered at Ao (call iV°), for simplicity, 
two other degenerate states, resonant at some wavelength 
Ao + A (or Ao — A), are formed. If these are coupled 
strongly enough, it can happen that a fortuitous mode 
splitting will push either the symmetric or the antisym- 
metric mode to fit into the band width of N°. In the case, 
when the original necklace at Ao is week, the new state 
can positively contribute to enhance the necklace trans- 
mission, thus increasing the order of N° by one. This is 
to say that the phase shift will count only one more res- 
onance, but still the intensity map will show extra nods, 
since the new resonance will appear with a double clump. 
We note, that in a random sample with resonant layers, 
such a situation cannot occur: the transmission spectrum 
in this case is always symmetric against c/Ao, therefore, 
the described contribution to the necklace N° will occur 
from both sides (Ao ± A) and consequently will result in 
a back repulsion of degenerate modes out of the necklace 
band. 

To conclude, we address through this study some pecu- 
liarities of the formation of optical necklaces, resonantly 
coupled degenerate modes, in one-dimensional random 
multilayer systems. We stress the importance of using 
the phase-jump method to reveal the exact order of neck- 
laces, in contrast to counting intensity or electric field 
clumps through the necklace. We show how the extra 
resonances are filtered out from the necklace when in a 
sample length L the number of degenerate in energy res- 
onances exceeds the optimal one. 
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